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Optomechanics concerns with the coupling between optical cavities and mechanical resonators. 
Most early works are concentrated in the physics of optomechanics in the small-displacement regime 
and consider one single optical cavity mode participating in the optomechanical coupling. In this 
paper, we focus on optomechanics in the extremely-large-amplitude regime in which a mechanical 
resonator is coupled with multiple optical cavity modes during the oscillation. We explicitly show 
that the mechanical resonator can present self-sustained oscillations in a novel way with limit cycles 
in the shape of sawtooth-edged ellipses and exhibit dynamical multistability. By analyzing the 
mechanical oscillation process and the accompanied variation of the optical cavity occupation, we 
develop an energy-balanced condition to ensure the stability of self-sustained oscillation. The effect 
of the mechanical nonlinearities on the dynamics of the mechanical resonator is also investigated. 

PACS numbers: 42.50.Wk, 05.45.-a 


I. INTRODUCTION 

Optomechanics has attracted much attention and is 
undergoing rapid development in the recent years in¬ 
is]. It concerns with the coupling between optical cavity 
modes and mechanical degrees of freedom via radiation 
pressure [5], optical gradient forces uni, photothermal 
forces m or the Doppler effect [H] . Various motivations 
have driven the development of this research field, such 
as detection of gravitational wave , more sensitive 

sensors of displacement, mass or force HMD, fundamen¬ 
tal studies of quantum mechanics |23H15| , preparation of 
macroscopic quantum state [26U2^ , quantum state trans¬ 
fer [5nil34| . novel nonlinear coupling quantum physics 
[55H57] . and quantum information processing [3Sj. Most 
of these motivations focus on the quantum level of op¬ 
tomechanical systems and require cooling the mechanical 
resonator as close as possible to its ground state, which 
has been widely investigated both theoretically [Mmol 
and experimentally |Ml - l48] . 

On the other hand, amplification of the mechanical os¬ 
cillation is very useful for both practical applications and 
fundamental researches, such as non-volatile mechanical 
memory [49) . synchronization of remote mechanical res¬ 
onators inn] and chaos dynamics [STjinn]- With a driving 
laser of frequency blue detuned with respect to the cav¬ 
ity resonance and power above a certain threshold, the 
mechanical resonator can run into self-sustained oscilla¬ 
tions mu mm [M]- An arbitrary tiny thermal fluctuation 
will be amplified into an oscillation with exponentially 
increasing amplitude, and finally be saturated into a sta¬ 
ble periodic oscillation. Self-sustained oscillation is of 
broad interests. It not only exists in optomechanical 
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system, but also in other systems such as a resonator 
driven by a superconducting single-electron transistor 
an ultracold atomic gas in an optical cavity m, 
or a metallic point contact deposited on a metallic spin- 
valve stack m- When optomechanical systems present 
self-sustained oscillations, dynamical multistability may 
emerge as theoretically predicted in |58j and experimen¬ 
tally explored in [59], which means that there may exist 
multiple possible stable oscillations at a set of fixed pa¬ 
rameters. 

Most early theoretical works on optomechanical self- 
sustained oscillation focus attention on the small- 
displacement regime. In this regime, it is sufficient to 
take only one single optical cavity mode into account par¬ 
ticipating in the coupling with a mechanical resonator, 
and the dependence of the resonance frequency of this 
mode on the displacement of the mechanical resonator 
can be treated linearly [5H]. The mechanical resonator 
conducts approximately sinusoidal oscillations at its in¬ 
trinsic frequency, and therefore, its limit cycles in the 
phase space are approximately elliptical. While in the 
large-amplitude regime, which can be realized in an op¬ 
tomechanical system driven by a high-power laser, the 
mechanical resonator can display different self-sustained 
oscillations with limit cycles mushroom-like in shape m- 

If the power of the driving laser is further increased, 
the amplitude of the mechanical oscillation can be com¬ 
parable with the wavelength of the laser, i.e., the system 
reaches the extremely-large-amplitude regime (ELAR). 
Multiple optical cavity modes of different orders may be 
excited and participate the coupling with the mechanical 
resonator during the oscillation. In this work, we will 
focus on optomechanics in the ELAR. We organize the 
paper in the following way. Sec. [^introduces the model 
and the Hamiltonian, and gives the dynamical equations 
of the system. In Sec. |III[ the self-sustained oscillation 
and dynamical multistability in the ELAR are studied by 
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analyzing the limit cycles of the mechanical resonator in 
the phase space. An energy-balanced condition is given 
to ensure the stability of self-sustained oscillation. The 
effect of the mechanical nonlinearities on the dynamics of 
the mechanical resonator is discussed in Sec. IV Finally, 
Sec. |V] gives a summary of this work. 


II. HAMILTONIAN AND DYNAMICAL 
EQUATIONS 


We consider a generic optomechanical system as shown 
in Fig. which is essentially a Fabry-Perot cavity with a 
fixed, partially reflecting mirror on one side and a mov¬ 
able, perfectly reflecting mirror on the other side. In 
this system, radiation pressure provides the dominant op¬ 
tomechanical coupling which is typically dispersive, im¬ 
plying that the primary effect of the moveable mirror is 
to shift the frequency of the optical cavity modes. The 



FIG. 1. (Color online) Schematic of a generic optomechan¬ 
ical system. An optical cavity of static equilibrium length 
I/O with a fixed, partially reflecting mirror on the left side 
and a moveable, perfectly reflecting mirror on the right side 
is driven by an external laser of frequency uji and power P. 
In the ELAR, multiple optical cavity modes of frequencies 
uj„ and decay rates can be excited and coupled with the 
mechanical resonator of intrinsic frequency ujm, mass m and 
damping rate 7 . 


cavity is driven by an external laser of frequency uji and 
power P. If the driving laser is turned off, the static equi¬ 
librium length of the cavity is Lg- The movable mirror in 
this model can be considered as a mechanical resonator 
of intrinsic frequency Wm, mass m and damping rate 7 . 
Usually, the displacement x of the mechanical resonator 
is assumed to be very small, so that it is sufficient to take 
only one single optical catity mode of frequency u}c{x) 
and decay rate k into account in the coupling with the 
mechanical resonator. In this small-displacement regime, 
iOc{x) is approximately equal to the first-order expansion 
around the static equilibrium position a; = 0 , 

uic.{x) « a;c(0) -I- uj'^{0)x = Wco - {uJco/Lq)x, 


where Wco = Wc(0). In the frame rotating at frequency 
w/, the Hamiltonian of the system can be written as 

H = h(uJco — 0Ji)a^a H- 1 — miJLx^ — hga^ax 

2m 2 

-t- (I) 

where a and are the bosonic annihilation and creation 
operators of the optical cavity mode, x and p are the po¬ 
sition and momentum operators of the mechanical res¬ 
onator, g = lOco/Lq is the optomechanical coupling, ul 
is the complex amplitude of the driving laser held which 
satishes = 2HiP/huJc. Here, without loss of general¬ 
ity, tti is set to be real. denotes the coupling between 
the optical cavity mode and the vacuum bath that leads 
to the decay rate n. refers to the interaction between 
the mechanical resonator and the thermal reservoir which 
is the cause of the damping rate 7 . 

If the optical cavity is driven by a high-power laser, the 
small-displacement assumption will be no longer valid 
and the system will reach the large-amplitude regime. 
The expansion of uJc(x) to the first-order will not be a 
good approximation. It is necessary to deal with uJcix) di¬ 
rectly without any approximate expansion. If the power 
P is high enough, the amplitude A of the mechanical os¬ 
cillation may be comparable with the wavelength A; of 
the driving laser. 


A/Xi ^ I, 


i.e., the system reaches the ELAR. In this case, not only 
ujc[x) should be dealt with directly, but also multiple 
optical cavity modes of different orders should be taken 
into account in the coupling with the mechanical res¬ 
onator during the oscillation. Concretely, each time the 
mechanical resonator passes through the positions that 
satisfy x Lq = n\i/2 {n S N), the nth-order optical 
cavity mode of frequency uJn(x) = WKclix P Lq) will be 
excited. So during a whole cycle of the mechanical os¬ 
cillation, a series of optical cavity modes will be excited. 
All these modes are coupled with the mechanical res¬ 
onator by radiation pressure. Here, we assume that the 
size of the mechanical resonator is much larger than the 
amplitude of the oscillation. In this situation, we can 
treat the mechanical resonator as a harmonic resonator 
and neglect the effect of mechanical nonlinearities (which 
will be considered in Sec. IV). Thus in the ELAR, the 


Hamiltonian of the system reads. 


H = h'^[uJnix) -UJl]a'ldn + ^ + 2 '^^mX^ 

n—1 


^(d„ -I- ^ -I- Hj, 


( 2 ) 


where a„ and ajj are the bosonic annihilation and cre¬ 
ation operators of the nth-order optical cavity mode of 
frequency ujn(x) and decay rates /v„. Heisenberg equa¬ 
tions of motion for operators d„, x and p can be easily 







3 


derived from the Hamiltonian above. In this paper, we 
aim to investigate the purely classical dynamics of the 
system, so we replace a„ by the complex coherent light 
amplitude a„(t), as well, x and p by their classical coun¬ 
terpart x{t) and p{t). We thus obtained the classical 
dynamical equations of the system as follows, 

/ TLTTC \ 

C^n ^ [ I Y I ^n^n-! (3) 

\x-\-Lo ) 

J_ oo 

2 li TITTC I |2 . /.s, 

For simplification, in the following, we assume all the 
are equal, namely, k„ = k. 


III. SELF-SUSTAINED OSCILLATION AND 
DYNAMICAL MULTISTABILITY 


Eqs. (|^ and Q are two coupled nonlinear differential 
equations, it is difficult to derive analytical solutions for 
them. Therefore we integrate these two equations numer¬ 
ically by using a fourth-order Runge-Kutta algorithm. 
We have considered the experimental feasibility [6TH64] 
and set the parameters as follows, the intrinsic frequency, 
mass, and damping rate of the mechanical resonator are 
iOm = 10^ Hz, m = 5 X 10“^® kg, and 7 = The 

wavelength and frequency of the external driving laser are 
A; = 1000 nm and uji = The static equilibrium 

length of the cavity is Lq = NXil2 and N — 10000. The 
decay rates of the cavity modes are Kn = k = lO^Wm, i.e., 
the system is in the unresolved sideband regime so that 
the optical cavity can respond quickly enough to the fast 
mechanical oscillation. It should be noticed here that, 
compared with Lg, the displacement of the moving mir¬ 
ror (i.e., the change of the cavity length) is usually very 
small even in the ELAR. Thus, we treat k as a constant 
rather than a function of x. 

Erom the numerical solutions of Eqs. <§-@. we plot 
the limit cycles in the phase space of the mechanical res¬ 
onator scanned by x and p = rax as shown in Eig. 
Each limit cycle corresponds to a stable self-sustained 
oscillation in the long-time limit. It is shown that the 
mechanical resonator can present self-sustained oscilla¬ 
tions in a novel way with limit cycles in the shape of 
sawtooth-edged ellipses. When P is relatively low, there 
is only one dynamical stable solution of the coupled Eqs. 
@-0 in the long-time limit, hence there is only one limit 
cycle in the phase space of the mechanical resonator as 
shown in Fig. Ha). The limit cycle is mushroom-like. 
Self-sustained oscillations with limit cycles of this shape 
have been theoretically studied m and experimentally 
observed [65j in optomechanical systems. 

When P is above 4.1W, there are multiple different 
limit cycles in the phase space of the mechanical res¬ 
onator as shown in Figs. Hb) -(h), meaning that the me¬ 
chanical resonator exhibits dynamical multistability. In 




x/Xi x/Xi 




x/Xi x/Xi 






FIG. 2. Limit cycles in the phase space of the mechanical 
resonator scanned by x and p = mx with different values of 
P. The parameters are: ui-m ~ 10^ Hz, m = 5 x kg, 

7 = Kn = K = lO^UJrn, Xl = 1000 Um, Lq = NXl/2 

and N = 10000. The values of P are: (a)P = 1 W, (b)P = 4.1 
W, (c)P = 7 W, (d)P = 10 W, (e)P = 11 W, (f)P = 12 W, 
(g)P = 15 W, (h)P = 17 W, as specified by the vertical black 
dashed lines in Fig. 


Fig.H we show the variation of the minimum and max¬ 
imum limit cycles respectively as P changes in the range 
from 4.1W to 18W. If the mechanical resonator initially 
reach a limit cycle in Fig. and P is then increased 
slowly, the mechanical resonator will oscillate with an ex¬ 
panding limit cycle changing shape in the way shown in 
Fig. Ha). While, if the mechanical resonator is initially 
trapped in a limit cycle in Fig. Hb), it will oscillate 
with an shrinking limit cycle as P decreased gradually as 
indicated in Fig. Hb). 

When the mechanical resonator moves forward (to the 
right as shown in Fig. H, the frequencies of the optical 
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FIG. 3. (Color online) The variation of the minimum(a) and 
maximum(b) limit cycles as the power of the external driving 
power changes in the range of 4.1W to 18W in which the 
mechanical resonator presents dynamical mnltistability. 


cavity modes oJn(x) = mrcj^x + Lq) = wKclix + N\il2) 
decrease. Conversely, when it moves backward, the 
frequencies increase. Every time the mechanical res¬ 
onator passes through the positions x = xu = kXi/2 
(k = ..., — 2 , — 1 , 0 , 1 , 2 ,...), the frequency of the optical 
cavity mode of order N + k meets the external driving 
laser frequency wi as shown in Fig. [^a). As a result, this 
optical cavity mode is excited and the photon number in 
the cavity becomes very large rapidly, accordingly, the 
radiation pressure exerted by light stored in the cavity 
increases rapidly. When the mechanical resonator goes 
away from these positions, the cavity occupation as well 
as the radiation pressure decreases rapidly to almost zero 
due to the large decay rates of the cavity modes. The 
variation of the photon number in the cavity accompa¬ 
nied with the oscillation of the mechanical resonator is 
shown in Fig. |^b). 

When the mechanical resonator is away from positions 
a; = Xfc, it approximately carries out damped harmonic 
oscillation, 

X = -xl,x - 'jx. ( 5 ) 

Every time the mechanical resonator passes through po¬ 
sitions X = Xk, radiation pressure in the cavity kicks it 
and dynamical equations of the system can be approxi¬ 





FIG. 4. (Golor online) (a)The dependence of the frequencies 
LJn of the optical cavity modes on the displacement x of the 
mechanical resonator. (b)The variation of the optical cavity 
occupation with the oscillation of the mechanical resonator. 
Here, the power of the external driving laser is P = 11 W. 
The other parameters are same as in Fig. 


mately written as. 


*5iv+fc {x a 

2 I ^9N+k I 
X = -Xrr,X-\ - 


- iaL - na 


N+k ’ 


( 6 ) 

(7) 


where = 47rc/[(A^ -I- k)X‘^] is the optomechanical 

coupling strength. It can be derived from Eq. Q that. 




ffLfc - Xk) + 


(a; - Xk) X 

(x - Xkf + 


( 8 ) 


Thus, the mechanical resonator satisfies the following 
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FIG. 5. (Color online) (a) Definition of average amplitude as 
A = + A‘^^^)l2. (b) Attractor diagram on a plane of 

A and P. Here, the parameters are same as in Fig. The red 
(dense) dots show the numerically exact results captured by 
solving the original coupled dynamical equations The 

black (sparse) dots are the approximately estimated results 
(see the Appendix) from the energy-balanced condition, Eq. 
( |12| |. The vertical black dashed lines specify some values of P 
distributed in different regions. Limit cycles for these values 
of P are plotted in Fig. 


equation of motion. 


X = —uj^x — ^x + 


hg 


2 

N + k 


{x — Xk)‘^ + 


jx - kX/2) X 
9 l+^ {x-Xkf 


(9) 


If the mechanical resonator moves forward, radiation 


pressure do positive work on it: 


VF+ 

N + k 




( 10 ) 


where x'^ is the velocity of the mechanical resonator just 
before it passes forward through the position x = Xk- As 
a result, the kicks lead to sharp accelerations of the me¬ 
chanical resonator. On the contrary, if the mechanical 
resonator moves backward, radiation pressure do nega¬ 
tive work on it. 


K SmK^^x'l ’ 


( 11 ) 


where x^ is the velocity of the mechanical resonator 
just before it passes backward through position x = Xk- 
So the mechanical resonator experiences sharp decelera¬ 
tions. As a reflection of these processes, there are some 
sawteeth on the limit cycles at positions x = Xk s& shown 
in Fig. 1^ If Xk is out of the oscillation range of the me¬ 
chanical resonator, we have Wt , = W~ ,, =0. When 
the total net work did by radiation pressure is balanced 
with the dissipative energy during one whole cycle: 


( 12 ) 


the mechanical resonator reaches stable self-sustained os¬ 
cillation. At some parameters, there may simultaneously 
exist multiple different possible stable oscillations that 
can satisfy this energy-balance condition. So the me¬ 
chanical resonator may exhibit dynamical multistability. 

To demonstrate dynamical multistability concisely, we 
define average amplitude as A = \/(^min + 
as shown in Fig. [§a), where Amin{Amax) is the mini- 
mum(maximum) amplitude of a limit cycle. It should be 
noticed that with the effect of radiation pressure the dy¬ 
namical equilibrium position x is shifted from the static 
equilibrium position a; = 0, but in the ELAR, the shift is 
small and can be neglected compared with the amplitude. 

Fig. [^b) shows the dependence of A on P, which can 
be considered as an attractor diagram. The red (dense) 
dots refer to the numerically exact results captured by 
solving the original coupled dynamical equations (©-Q- 
The black (sparse) dots denote the approximately esti¬ 
mated results (for details, see the Appendix) from the 
energy-balanced condition, Eq. The discrete aver¬ 

age amplitudes of the mechanical resonator reveals that 
the energy-balance condition leads to an amplitude lock¬ 
ing effect. When P is above 4.1W, the mechanical res¬ 
onator exhibits dynamical multistability, so A can take 
multiple values for a fixed P as shown in Fig. [^b). 


IV. MECHANICAL NONLINEARITIES 

In ELAR, if the size of the mechanical resonator is not 
very large, the stress induced in it during the oscillation 
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FIG. 6. (Color online) The dependence of the average ampli¬ 
tude A on a. Here, the power of the external driving laser is 
set to be P = 7 W. The other parameters are same as in Fig. 

H 


may be strong and the resulting mechanical nonlinearities 
[55H55] may be non-neglectable. In this situation, the dy¬ 
namics of the mechanical resonator can be described by 
the following Duffing equation, including optomechanical 
coupling, 


mrc 


X = -w^(l + ax^)x + — y] , , r ^2 

n—l ^ ' 


Ql-y] 


-IX, 


(13) 


where a is the cubic nonlinear constant and has the units 
of m“^. If the spring constant of the mechanical res¬ 
onator weakens with increasing amplitude, the value of 
a is negative. On the contrary, a positive value of a in¬ 
dicates that the spring constant stiffens with increasing 
amplitude. Fig. (§ shows the dependence of the average 
amplitude A on a. When a < 0, the average amplitude 
A changes gently with a, meaning that for a weakening 
mechanical spring, the mechanical nonlinearities have a 
weak effect on the mechanical oscillation. For a stiffening 
mechanical spring (a > 0), the effect is also very weak 
when a is small. Only when a is large, the effect becomes 
strong. As a increases, A increases sharply and the dis¬ 
tribution of the dynamical multistability is affected. 


V. SUMMARY 

We have studied classical dynamics of a generic op¬ 
tomechanical system in the ELAR in which a mechan¬ 
ical resonator is coupled with multiple optical cavity 
modes. It has been shown that the mechanical resonator 
can present self-sustained oscillations in a wide range of 
parametric space. Instead of sinusoidal oscillations in 


the small-displacement regime, the mechanical resonator 
display more complicated oscillations in the ELAR with 
limit cycles in the shape of sawtooth-edged ellipses. We 
have analyzed the process of the mechanical oscillation 
and the accompanied variation of the optical cavity occu¬ 
pation. Based on these, we derived an energy-balanced 
condition of stable self-sustained oscillation. We have 
demonstrated that the mechanical resonator may ex¬ 
hibit dynamical multistability, which can be explained 
by the fact that there may simultaneously exist multiple 
different stable oscillations that can satisfy the energy- 
balanced condition. The effect of the mechanical nonlin¬ 
earities on the dynamics of the mechanical resonator has 
also been discussed, and it has been shown that the ef¬ 
fect is weak in a wide range of values of the nonlinearity 
parameter a. Only when a is positive and very large, the 
effect becomes strong and may has an influence on the 
distribution of the dynamical multistability. 

Dynamical multistability may find applications in sen¬ 
sitive force or displacement detections and memory stor¬ 
age. In the parametric space of our calculation in this 
paper, the required power of the external driving laser 
is very high. However, with smaller mass, intrinsic fre¬ 
quency or damping rate, the mechanical resonator can 
reach the ELAR driven by a much weaker laser field. 
It is possible to observe self-sustained oscillation demon¬ 
strated in this paper within the reach of current experi¬ 
mental technology. 
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APPENDIX 


The energy-balanced condition, Eq. (12) can be rewrit¬ 
ten as an equation of two unknown quantities Amin and 
Amax as follows: 


(w+ +w~ \ 

/ ^ y N + k ' N + k J 


{ — Amin<Xk<A.max) 


= E^{A 


7 yx^mini -^maxJi 


(14) 


where and are expressed by Eq. (10) and 

0 - It should be noticed that Amin and Amax are not in¬ 
dependent of each other, but are correlated. Concretely, 
if the mechanical resonator moves forward with initial po¬ 
sition and velocity {x,x) = (—Amm,0), it approximately 
experiences the following motion processes. It carries out 
damped harmonic motion obeying Eq. <§ until it reach a 
state {xk,x'^)- Next, it experiences a sharp acceleration 
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rapidly and its state undergoes an abrupt change: 


{xk,ik) 



+ 3h'^0‘t9N+k'n- 
V A; / rriK 4m^ 


(15) 


which can be derived from Eq. (10). Then the mechan¬ 
ical resonator alternately experiences damped harmonic 
motions and sharp accelerations at the positions x = Xk 
until it reach a unique and determinate state {A^axfi)- 
Similarly, when the mechanical resonator moves back¬ 


ward, it alternately experiences damped harmonic mo¬ 
tions and the following sharp decelerations at the posi¬ 
tions Xk' 


{Xk,X+) 



V fc / rriK. 


(16) 


By numerically simulating these motion processes, we can 
get values of VEjv+fc+, ^N+k testing whether 

Eq. (13) is satisfied, we can get Amin, A^ax, and then 
A of stable self-sustained oscillation. 
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